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' O ' Abstract. - We show that a generalized charge SU{2) symmetry of the one-dimensional (ID) 

C , Hubbard model in an infinitesimal flux generates half-filling states from metallic states which 

lead to a finite charge stiffness D(T) at finite temperature T, whose T dependence we study. Our 
results are of general nature for many integrable quantum lattice systems, reveal the microscopic 
mechanisms behind their exotic T > transport properties and the interplay with T = Q quantum- 
phase transitions, and contribute to the further understanding of the transport of charge in systems 
^ ' of interacting ultracold fermionic atoms in ID optical lattices, quasi-lD compounds, and ID nanos- 

. tructures. 

d 

I There has been a recent increasing interest in the unusual transport and spectral prop- 

. erties of systems of interacting ultracold fermionic atoms in one-dimensional (ID) optical 

lattices [1] and a renewed interest in those of quasi-lD carbon nanotubes, ballistic wires, and 
quasi-lD compounds [2,3]. Quantum effects are strongest at low dimensionality leading to 
' unusual phenomena such as charge-spin separation at all energies [3] and persistent currents 

^ . in mesoscopic rings [4]. The ultracold atom technology can realize a ID closed optical lattice 

experimentally with a tunable boundary phase twist capable of inducing atomic "persistent 
currents" [1]. Thus, the further understanding of the microscopic mechanisms behind the 
transport of charge in low-dimensional correlated systems and materials is a topic of high 
scientific and technological interest. Recently there has been an enormous grow in the lit- 
erature on the quantum phase transitions in such systems [5]. For instance, the relation 
of the unusual T > charge-transport properties observed in low-dimensional correlated 
systems to their T — quantum phase transitions is also a problem which is far of being 
fully understood. 

For the ID Hubbard model [6] or any other interacting electronic quantum system defined 
in a ID lattice of Na sites of length L = aNa = Na {ks ~ h = a = — e = 1) the T > 
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ideal conducting and insulating behaviors are defined in terms of the real-part of the optical 
conductivity for finite temperatures T > 0, which reads ct{u>) = 2-k D 5{lo) + Uregi^)- Here 
the charge stiffness is given hy D = D{N,T) = d/Z where Z = Z{N,T) = J^mPm is the 
partition function of the A^-electron system, the m summation runs over the eigenstates 
|m) of energy Em, and Pm = e"^""/-^ are the corresponding Boltzmann weights. At T = 
the charge stiffness is related to the boundary energy of the lattice model under twisted 
boundary conditions [7] and for T > d = d{N, T) can be related to the thermal average of 
curvatures of energy levels subject to an infinitesimal flux (j) (in units of 27r/L and L oo), 

'^=L22Pm-^^^ - (1) 

At T > a normal conductor is such that D{N, T) = and cr(a; ^ 0) > 0, whereas an ideal 
conductor is characterized by D{N,T) > 0. Furthermore, an insulator can develop into 
a normal conductor, become an ideal conductor, or display ideal insulating behavior such 
that -D(T) = and a{iv 0) = 0. ID integrable systems whose Hamiltonians commute 
with the current operator behave as ideal conductors and ideal insulators for the densities 
corresponding to the T = metallic and insulating quantum phases, respectively. There 
is strong numerical evidence that the same occurs concerning the T > ideal conducting 
behavior for densities referring to T = metallic phases of ID integrable quantum systems 
whose Hamiltonian does not commute with that operator, such as the ID Hubbard model 
[8,9]. However, whether for the densities corresponding to the T = insulating quantum 
phases such systems behave as ideal insulators remains an open question. In contrast, 
the non-integrable ID interacting systems are generic conductors and activated ones in the 
metallic and insulating phases, respectively. 

In this Letter we clarify such an issue by showing that a generalized charge SU{2) 
symmetry of the ID Hubbard model in an infinitesimal flux (f) generates half-filling states 
from metallic states which lead to a finite half-filling charge stiffness D = D{Ng,,T) at 
finite temperature T. We study its T dependence by characterizing its different behaviors 
in three well-defined temperature regimens. Interestingly, the above half-filling states do 
not belong to the integrability Hilbert subspace where the model is solvable by the Bethe 
ansatz [6]. The latter subspace is spanned by the lowest- weight states (LWSs) of both 
the r]-spin and spin SU{2) algebras [10]. When defined in such a subspace, the quantum 
problem has an infinite number of conservation laws [11]; as in the corresponding classical 
systems, the occurrence of an infinite number of such laws is behind the integrability of 
the quantum models. Moreover, we present numerical evidence that when acting onto the 
integrability subspace the ID Hubbard model is an ideal insulator. Indeed, we find that 
in the thermodynamic limit the half-filling states that span such a subspace do not carry 
charge current, consistently with the general predictions of Ref . [8]. Importantly, our results 
provide useful information about the microscopic mechanisms behind the interplay of the 
T = Mott-Hubbard insulator - metal quantum-phase transition with the T > thermal 
properties in ID correlated systems and new insights into the exotic microscopic processes of 
charge transport in systems of interacting ultracold fermionic atoms in ID optical lattices, 
carbon nanotubes, ballistic wires, and quasi-lD compounds [1-3]. 

The ID Hubbard model in an infinitesimal flux (p reads, 

H{<f>) = -tY,W^ 4, a^i+i, <^ + h.c] +UD., (2) 

where D- = J2j{{fij,-t - 1/2] ['^j,i - 1/2] + 1/4}, c]_^ and Cj^a (and below and Cfe^o-) 
are spin-projection a =t, i electron operators at site j = l,2,...,A^a (of momentum /c), 
the number Na is even and very large, and fij^^ = Cj^^Cj^c^. We use periodic boundary 
conditions and consider zero spin density [n-f — nj =0 and electronic densities n = [n-f+ni] 
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and corresponding hole concentrations 6 = 1 — n in the range < S < 1, where = N„/L = 
Na/Na and N„ such that N = [N^ + A/'J is the number of a electrons. The following ^ = 
properties of the Hamiltonian (2) play an important role in our study: (i) Its zcro-cncrgy 
level refers to the 6 = insulating ground state, Eo{Na) = 0, and is located in the middle 
of the Mott-Hubbard gap [6], Amh = [i-itf/U] du:\\/J^ - l/sinh(27rta;/C/)], such that 
Amh w [8\/[7t/7r]e-2'^*/'^ for V jt « 1 and Amh ^ U - At for U/t » 1; (ii) The ground- 
state energies associated with even values of N < Na are such that Eo{N—2) > Eo{N) where 
£;o(-^V) changes from [47V„V7r][l-cos(7r^/2)] for U/t«l to [SNaU/2]{l-[Usm{Tr5)/Tr5U]} 
for U/t >> 1; (iii) The T = chemical potential 2/x(iV, 0) is an increasing function of S, 
changing from 4tsin(7r5/2) for U/t « 1 to C/ — 4tcos(7r5) for U/t » 1 and is given 
by 2/i(A/',0) w Amh + S^n'^/ml and 2/i(0,0) = U + At > Amh for 5 small and 6 = 1, 
respectively, and all values of U/t; (iv) is a charge mass given in Eq. (48) of Ref. [12] 
whose limiting values are ^ as U/t — *■ and l/2t as U/t oo. 

For (j) = the Hamiltonian (2) commutes with the three generators of the ry-spin SU (2) al- 
gebra [10]. While for (/> > it docs not commute with the two off-diagonal generators of that 
algebra, a result which plays a central role in our study is that for (j) > it commutes with the 
three generators % = -^ Ej [1 - ^j,T -nj,i]^'nl = J2k 4-0-^, i 4, t = e^^'^'^"'^^ ^ i 4 T ' 
and fj^ = Cfc^ I CT^-4,-k, [ = Y^j ('^''^^'"^'^''cj^ ^ Cj^i of a generalized charge SU{2) algebra. 
Since for > the number of electrons N and r^-spin value r] remain good quantum numbers, 
the diagonal generator i)^ is independent of (j). Furthermore, note that as ^ — > the above 
off-diagonal generators of the generaUzed algebra become those of the ry-spin algebra [10]. 
For (j) > the Hamiltonian (2) remains integrable. However, as for ^ = 0, the generalized 
Bethe-ansatz equations given in Eq. (10)-(12) of Ref. [9] for (/> > refer to the Hilbert sub- 
space spanned by the LWSs of both the ?7-spin and spin algebras such that 2ri = [Na — N] 
and 25* = [iV-f — iVj , where S is the spin. Thus, for half filling such an integrability subspace 
is spanned by the ?? = 0, /y-spin singlet states. The full Hilbert space contains another sub- 
space which is spanned by energy eigenstates generated from application onto the LWSs of 
the off-diagonal generators of the generalized algebras. Since here we are mostly interested 
in the transport of charge, for simplicity we limit our study to the charge and ?7-spin degrees 
of freedom, yet it is straightforward to extend our analysis to the spin degrees of freedom. 
Let [to, N, (j)) be a A''-electron energy eigenstate of the Hamiltonian (2) for > 0. such that 
2r] > [Na — N] . This state does not belong to the subspace where the model is integrable but 
can be generated from the corresponding A'''-electron LWS such that A''' = [A^^ — 2?7] < A'' 
as follows, 

|m',Ar',<^), (3) 
where the normalization constant reads C = 6^^ i [;Va-Ar] + nl'=i^'^''~^^^[2r/+l-^]- Since the 



|m, Ar,( 



i 

77^ 
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Fig. 2: The degenerated energy spectrum and charge current spectrum in units of t of the half- 
filling ry-spin-triplet states (left current spectrum) and corresponding singlet states (right current 
spectrum) for U/4t = 10 and iV = iV„ = 66. 



Hamiltonian (2) commutes with the generators fjjp and 77^, the energies Em{(f>) and Em' {4') 
of the above states |m, N, cp) and |m', N', (p) are identical. 

The microscopic mechanism described in the following implies that, in spite of the model 
integrability, for finite values of U/t the ^ = charge stiffness is finite for finite T values 
and vanishes both as T ^ and T ^ 00. For each energy eigenstate |m, N, cj)) of electronic- 
number N and 2rj > [Na — N] there is one and only one Ty-spin LWS \m',N',(j)) with 
the same r] value, N' < N, and 2t] = [Na — N'] whose energy is such that Em'{4>) = 
Em{4>), where Em{(l)) is that of \m,N,(j)). In turn, there are 2r] states \m,N,(j)) such that 
N = Na — 2ri + 2j with j = 1, 2, 2?7 corresponding to the same r]-spin LWS \m' , N' , (p) . 
Thus, Z{N,T) = 6n,n^ + Zlws{N,T) + Zr{N,T) = dN,N^ + En'=o Zlws{N' ,T) (and 
d{N,T) = dLws{N,T)+dR{N,T) = E^'^o^LWsiN' ,T)), wheve 5 n,n^+Zlws {N,T) (and 
dhwsiN, T)) corresponds to the contributions from all the rj-sp'm LWSs of electronic number 
N and Zr{N, T) = ^^.=0 Zlws{N', T) (and dR{N, T) = Y^^V^q dLws{N', T)) is such that 
the A''' summation refers to all LWSs with even electronic numbers N' = 0,2,4, ...,N — 2 
from which one can generate all non-LWSs with A'^ electrons. (The above contribution S^.Na 
refers to the insulating ground state.) Such results imply that ma,x d{N,T) = d{Na,T) and 
thus if d{N,T) > for iV even, AT < A/-„, and T > then d{N,T) > ior N = Na and 
r > 0. 

Without the construction of the states (3) for </> > the above results could not be 
reached. Their use reveals that a parameter which plays key role in our study is, 

a{AN,T) = Z{AN,T)/Z{Na-AN,T); Z{AN,T) = ^ Zlws{N',T), (4) 

N'=Na-AN-2 

where Z{Na — AN,T) is the partition function for N = Na — AN and the N' sum runs 

over even values only. For finite values of T and U/t one has that a{AN,T) < 1 for 
AN < ANo = ANo{T), where the critical value AiVo(r) = 2,4,6,... increases with T, 
being such that ANo{T) ^ 2 as T ^ and ANo{T) ^ Na as T ^ 00. Moreover, 
min{a(A7V, T)} = a{2, T), max{a(2, T)} = a{2, 0) = 1, and a{2, T) < 1 for T > 0. 

Three typical regimens correspond to (a) 2iJ,{Na,T)/T 00 when T « 2/u(iVa,0) = 
Amh, (b) 2ij,{Na,T)/T Ki for intermediate temperatures T 2fi{Na,0) = Amh, and (c) 
2/i(0,T)/T ^ -00 when T ^ 00. We find from the above results that for U/t > the 
half-filling charge stiffness is of the general form, 

-<"--) = ??^-'-«.-'. 

where F{AN,T) = a{AN,T) + l/Z{Na - AN,T) and 

7(AAr, T) = [d{Na, T) - d{Na - AN, T)]/Z{Na, T) , (6) 
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Fig. 3: The charge current spectra in units of t of the half-filling 77-spin-singlet states of Figs, f and 
2 for U/4t = 1 and N = Na = 18, N = Na = 30, and N = Na = 46. 



is such that j{AN,T) > for all AN 2,4, ...,Na. For the regimen (a), e'^"^"/'^ » e"°° 
and thus also e-2/'(^'.o)/J' « 6"°°, because 2fi{N',0) > Amh for < iV' < Na. Thus, one 
has both that , 

ZLWsiN\T) _ 2M(iV',0)/T _ -00 . dLWs{N',T) _ 2,,(iV',0)/T _ -00 

^ws(A^'-2,T) ~ ~ rfws(A^'-2,T) ~ 

(7) 

and only the term dLws{Na - 2,T) contributes to dR{Na,T) = J2n'^=o dLWs{N' ,T). We 
then find from use of Eq. (5) that, 

D{Na,T)(xe-^""/^^T/ml; T « Amh ■ (8) 

In turn, as the temperature increases and reaches the regimen (b) such that T « Amh , one 
finds that e2''(^''^)/^ « 1 for an increasing number of Zlws{N', T) and dLws{N', T) terms 
of the expressions ZR{Na, T) = J^^^J Zlws{N', T) and dniNa, T) = ^wCo dLws{N', T),| 
respectively, so that ZR{Na,T) and dR{Na,T) have significant contributions from a maxi- 
mum number of such terms. Moreover, So = A7Vo(T) /Na is now finite and e2A'(Wa-AAfo,T)/T ^| 
1 implies that Z{Na-ANo, T) » 1 and thus F(A7Vo, T) < 1 or F{ANo, T) « f . Therefore, 
since for 5q finite the system is an ideal conductor for TV = Na — ANq [8,9], D{Na,T), Eq. 
(5), is finite for the regimen (b), Furthermore, the achieved maximum contribution from 
the whole set of above terms leads to a maximum value of D{Na, T) for T « Amh- For in- 
creasing values of T, D{Na, T) becomes a decreasing function of T and finally when T ^ 00 
the regimen (c) such that the two ratios of Eq. (7) reach the values e~'^^^^ >t)/t ^ gOo 
g2n{N ,T)/T goo^ respectively. We then find from use of Eq. (5) that, 

D{Na,T)^l/T; T^oo, (9) 

which is the expected behavior of systems on a lattice. The overall charge-stiffness T de- 
pendence is thus similar to that of the 6 = dimerized model plotted in Fig. 6 of Ref. [9] , 
with the dimerized gap A replaced by Amh- The temperature T rs Amh associated with 
the maximum value of D{Na,T) is an increasing function of U/t. For U/t ^ 00 one has 
that T « Amh ~ C/ — 4<: — > 00 and thus D{Na,T) is given by Eq. (8) for all ranges of 
finite T values. Thus, for U/t — > 00 the regimens (b) and (c) are never reached for fi- 
nite T values. In the opposite limit, when upon decreasing the value of U/t one reaches 
Tk. Amh^ [8VC^/7r]e-2'^*/^, the temperature range above T = where D{Na, T) IS given 
by Eq. (8) becomes smaller and smaller, until it shrinks as both Amh and — > 
and the limit U/t = is reached. In that limit a maximum value D{Na, 0) = 2f/7r occurs for 
T = Amh and thus D{Na,T) becomes finite and largest at T = 0. Thus, the regimen 
(a) is absent for U/t — > 0, which implies a metallic ground state for U/t = and S = 
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Fig. 4: The square root of the summation of the currents of the half-filhng T^-spin-singlet energy 
eigenstates with one charge string of length two (left) and two strings of length one (right) divided 
by the number of these states for U/4t = 10 as a function of 1/L. The dashed lines correspond to 
the expected behavior for 1/L — > 0. 



and ideal conducting behavior associated with the regimen (b) for finite temperatures. Our 

above results then reveal the microscopic processes behind the T = half-filling U/t ^ 
Mott- Hubbard insulator - metal phase transition. Indeed, for finite U/t values D{Na,T) 
displays a maximum value at T Amh which corresponds to the merging of the regimen 
(b) and vanishes both for T ^ and T ^ oo as given in Eqs. (8) and (9), respectively. 
However, for U/t ^ the maximum value is shifted to T = as Amh 0, which is 
behind the T = quantum phase transition. The above mechanism reveals that the finite 
values D{Na,T) > for T > and U/t > result from half-filling states generated by the 
generalized charge SU (2) algebra from metallic states. 

For L — > oo the charge stiff'ness can be expressed as D = [1/2TZL] J2mPm (im)^ in terms 
of the expectation values jm = {m\j\m) = —[dEm{4>) /d4>\\0^ q of the current operator j and 
thus D = D{N,T) vanishes if jm = for all energy eigenstates [9] . In order to confirm 
that the ^ = finite charge currents jm are carried by states generated from metallic states, 
wc provide numerical evidence that the (5 = ry-spin singlet energy eigenstates which span 
the integrability subspace carry no charge current and thus that dLwsi^aiT) = 0. Our 
numerical results also confirm that the ^ = 0, 77 > states carry finite current, in agreement 
with our above analysis. Wc have used the L » 1 Bcthe-ansatz equations of Rcf. [9] with 
L finite and checked that our finite-size results are in excellent numerical agreement with 
all known quantities. In Figs. 1 and 2 we plot the ^ = energy spectrum as a function 
of the momentum of the rj = I, 77-spin-triplet energy eigenstates and of the corresponding 
77 = 0, ry-spin-singlet states whose energy spectrum is degenerated with that of the former 
states for two values of U/At. Such 77 = states have one charge string of length one [9]. 
These figures also display the charge-current spectrum of these two types of states. Note 
that the current value is finite for the //-spin triplet states, in agreement with our above 
analysis, whereas for the 77 = states it is very small for the finite system and vanishes as 
L ^ 00. In figure 3 we show the current spectrum of the same 77 = states for U/At = 1 
and different values of Na- The data of the figure confirm that the current values decrease 
for increasing values of L and our extrapolation for L — > 00 reveals that they vanish in that 
limit. Finally, we study the current expectation values of all (5 = 0, 77 = 2 states and of the 
corresponding 7/-spin-singlet states whose energy spectrum is degenerated with that of these 
states. There are two types of such 77 = states, which have one charge string of length two 
and two charge strings of length one. respectively [9]. We consider the square root of the 
summation of the currents of these two classes of states. The results are plotted in Fig. 4 for 
U/At = 10 as a function of 1/L and indicate that the currents of these states also vanish as 
l/L ^ 0. A similar analysis for other values of U/4t and random samples of 77-spin-singlet 
states whose energy is degenerated with that of (5 = 0, 77 > 2 states also led to vanishing 
currents as L — > 00. 
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In this Letter we have shown that a± N = Na the ID Hubbard model is not an ideal 
insulator because a generalized charge SU{2) symmetry of the model in an infinitesimal 
flux generates half-filling states from metallic states that carry finite current. Moreover, we 
found evidence that when acting onto the Hilbert subspace spanned by the energy eigenstates 
associated with integrability the ID Hubbard model is at half filling an ideal insulator. Our 
results also reveal the microscopic mechanisms behind the interplay of the half-filling T = 
Mott-Hubbard insulator - metal quantum-phase transition with the model unusual T > 
properties and contribute to the further understanding of the transport of charge in systems 
of interacting ultracold fermionic atoms in ID optical lattices, quasi-lD compounds, and ID 
nanostructures. 
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